Dynamics of scaled norms of vorticity for the three-dimensional 
Navier-Stokes and Euler equations 

J. D. GibboiQ 

Department of Mathematics, Imperial College London SW7 2AZ, UK 

Abstract 

A series of numerical experiments is suggested for the three-dimensional Navier-Stokes and Euler equations 
on a periodic domain based on a set of L 2m -norms of vorticity Q. m for m > 1. These are scaled to form the dimen- 
sionless sequence D m = (GJq 1 Q. m ) a "' where CJo is a constant frequency and a„, = 2m /(Am — 3). A numerically 
testable Navier-Stokes regularity criterion comes from comparing the relative magnitudes of D,„ and D m+ \ while 
another is furnished by imposing a critical lower bound on J ' D m dr. The behaviour of the D m is also important 
in the Euler case in suggesting a method by which possible singular behaviour might also be tested. 

To appear in the Procedia IUTAM volume of papers Topological Fluid Dynamics II. 



1 Introduction 



The challenges that face those concerned with the numerical integration of the three-dimensional incompressible 
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Euler and Navier-Stokes equations for a velocity field u(x, t) on a 3D periodic domain V = [0, L] 3 



Du 

= vAu-vp divu = (1) 

Dt F 

(v = for the Euler equations) are also reflected in the challenges faced by analysts in their attempts to understand 
the regularity of these equations. The best known and most effective result in which analysis has guided numerics 
is the Beale-Kato-Majda (BKM) theorem [QIU which says that solutions of the three-dimensional incompressible 
Euler equations are controlled from above by J ' \\co\\aadT. If this is finite then no blow-up can occur at time t. 
Moreover, any numerical singularity in the vorticity field of the type ||ft)||oo ~ (to — t) ~ p must have p > 1 for the 
singularity not to be a numerical artefact. The BKM criterion has become a standard feature in Euler computations : 
see the papers in the special volume |3J]. 

This present paper is concerned with regularity criteria that form a consistent framework for the Euler and 
Navier-Stokes equations and which are testable numerically. In both cases it is often not clear when a numerically 
observed spike in the vorticity or strain fields remains finite or is a manifestation of a singularity. It is well known 
that monitoring the global enstrophy, or //i-norm ||ffl||2, pointwise in time determines Navier-Stokes regularity 
Hi , while the monitoring of ||fl)||oo likewise determines the fate of Euler solutions. However, the range of V- 
norms between these may be useful. The basic objects are a set of frequencies based on L 2m -norms of the vorticity 
field co = curl u 

, \ l /2m 

&m(t) — [ L~ 3 / \co\ 2m dV) l<m<°o. (2) 



Holder's inequality insists that Q.„, < £l m +i- The Navier-Stokes and Euler equations are invariant under the re- 
scaling x' = Ex ; t ' = e 2 t ; u = eu' ; p = e 2 p'. If the domain length L is also re-scaled as L' = eL then £2 m re-scales 
as i2„, = £ 2 £l' m as expected. If, however, L is not re-scaled but kept fixed then Q. m re-scales as 

£!*■ = eQ> , (3) 

where 

""I = ~A 7 • W 

Am — 3 

It turns out that this strange scaling is particularly important and provides a motivation for the definition of the set 
of dimensionless quantities 

D m (0 = (d 1 £2 m ) a '", l<m<oo (5) 
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where a,\ = 2 and = 1 /2. For the Navier-Stokes equations the frequency GSq is easily defined as G7o = vL~ 2 . 
The case of the Euler equations is more difficult as there is no obvious material constant to replace v in the 
definition of flJo- A circulation F has the same dimensions as that of v but it must be taken around some chosen 
initial data: for instance, in fl, initial data was taken to be a pair of an ti -parallel vortex tubes, in which case T 
could be chosen as the circulation around one of these. For a discussion of the variety of conclusions that can be 
drawn from numerical experiments see fi0.ll.SE Ill- 



No proof exists, as yet, of the existence and uniqueness of solutions of either the 3D Navier-Stokes or Euler 
equations for arbitrarily long times. A time-honoured approach has been to look for minimal assumptions that 
achieve this result 14j,|5[]. In fact, much of what is known about solutions of both the 3D Navier-Stokes and Euler 
equations is encapsulated in the sequence of time integrals 



Jo 



D 2 m dx (6) 

o 



based on the continuum lying between 

£>i = flJ 2 ZT 3 ^ \co\ 2 dV ... — » ... Doo= (gJq'IMIoo) 172 . (7) 

A well-known time-integral regularity condition for the Navier-Stokes equatons is that the first in the sequence in 
© should be finite Q, S) : that is D 2 dx < °°. In contrast, the boundedness of the last in the sequence in (O 
at m = °° is exactly the Beale-Kato-Majda criterion J WcoW^dx < °° for the regularity of solutions of the Euler 
equations HUH. 

For three-dimensional Navier-Stokes turbulence, it has to be admitted that arbitrarily imposed regularity as- 
sumptions, such as f^D^dx < °°, while mathematically interesting, have little foundation in physics: see the 



discussion of this point in 11211 . However, what is known, without any assumptions, is that weak solutions (in the 



sense of Leray lfl3IO obey the time integral 1 14] 



D,„dx<c(tRe i + r] l ) , (8) 



where 771 is a constant depending upon D m (Q). This result plays two roles. In $3] it is shown that it leads to 
a definition of a continuum of inverse length scales LA" 1 the upper bound on the first of which is the well- 
known Kolmogorov scale proportional to Re 3 / 4 . The more general upper bound is discussed in fj3]and is given by 
LA" 1 < cRe i l 2a "'. Thus the A,„ for m > 1 correspond to deeper length scales associated with the higher L 2 "'-norms 
of vorticity implicit within the D m . 

In addition, the magnitude of the bounded time integral in ([8]) is also significant. Let us consider whether the 
saturation, or near saturation, of this time integral plays any role in the regularity question. In the forced case it has 
been shown in lfl2ll that if a critical lower bound is imposed on this time integral in terms of the Grashof number 
Gr then this leads to exponential collapse in the D m (t). In fact boundedness from above of any one of the D m also 
implies the boundedness of D\ which immediately leads to the existence and uniqueness of solutions. While it 
can be argued that the imposition of this lower bound is physically artificial the result is nevertheless intriguing 
because it suggests that if the value of the integral ([8]) is sufficiently large then solutions are under control, which 
is surely counter-intuitive. Once it dips below this critical value then regularity could potentially break down. 
The importance of this mechanism lies in the role it may play in understanding the phenomenon of Navier-Stokes 
intermittency. This is discussed in |J5] where the critical lower bound is expressed in terms of the more physical 
Reynolds number Re 



JS, 



[t Re i0 "' + T] 2 )< J D m dx, < 5m < 1 . (9) 

The range of 8,„ is estimated and it is shown that 8 m \ \ for large m, thus allowing considerable slack between the 
critical lower bound in © and the upper bound in ([8]). 

The D m are comparatively easy quantities to calculate from a numerical scheme and is thus it is worth exploring 
whether regularity criteria can be gleaned from the relative magnitudes of the D„, or their time integrals. This is 
the task of this paper. In Sj2] the Euler equations are discussed in these terms where two versions of a numerical 
experiment are suggested for testing singular or non-singular behaviour. 
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2 The incompressible 3D Euler equations 

Whether the three-dimensional Euler equations develop a singularity in a finite time still remains an open problem 
but a variety of super-weak solutions have recently been shown to exist 

EE [HE! EH!. Let r be the 
circulation around some chosen initial data such that GJo = IX , as discussed in SjT] This defines GJq within the 
definition of D m . 

Proposition 1 Provided solutions of the three-dimensional Euler equations exist, for 1 < m < °° the D m formally 
satisfy the following differential inequality 

I D \ 

b m <c m a { -^pj D 3 m , § m = !(4m+l). (10) 

Proof : The time derivative of the Q. m obeys 

2mL 3 n2'"- 1 n m = — f \co\ 2m dV 
dt J 

n [ (0 2m \Vu\dV 



< 2m \co\ 2m dVj \(0\ 2{ - m+ ^dVy m+1 yj \Vu\ 2{ - m+V > dV 

< 2mL 3 Cl , m n™+|£2™ (11) 



where we have used |V«|| p < c p ||o)||p, for 1 < p < °°. Thus it transpires that 



n m <c Um ^-^j Q. z m , l<m<°°. (12) 

Note that the case m — °° is excluded : it was shown in [1] that a logarithmic H3 = Jy \ V^u^dV factor is needed 
such that ||Vw|| 00 <c]|fi)||oo(l +lnH 3 ). 

We wish to convert inequality (V2\ to one in terms of D„, 

n \ m+l /n \ ~£r~ ( m+1 s+lU 1 / n \ 4( 4m +!) 

£i m = vJ^Y D * +l am) am =m ^-Y D 2 n 



having used the fact that 

-^ = 2, A„ = jm(m+l). (13) 

, Ct;n+ 1 G-m / 

Substitution into (fT2l completes the proof. ■ 
There are now at least two interesting routes for the integration of fllOt . 

1. Firstly if a finite time singularity is suspected then divide (fTOb by £> 3 „~ £ with < £ < 2 to obtain 



where 



/' / D 1 \ 



For instance, for e = we have 



(15) 



* {t) * P.M-*-FM (16) 
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where 

F lfi {t)=2c m m J{-^-J dz. (17) 

A singularity in the upper bound of inequality ( TTol ) is not necessarily significant. What is more significant is 
whether the solution tracks this singular upper bound. This suggests the following numerical test : 

(a) Is Fifi linear in t ? 

(b) If so, then test whether 

D 2 m {T c , m -t)^C m with T c ,m^T c (18) 

uniformly in m. If such behaviour occurs it suggests, but does not prove, that the D,„ may be blowing 
up close to the upper bound. 

2. If exponential or super-exponential growth is suspected then divide ( TTOb only by D m (the case £ = 2) and 
integrate 



D m <D m (/o)expGJ / ( ) ' D l d ?- (W) 



The rate of growth of the integral with respect to t is of interest. Does it remain finite for as long as the 
integrated solution remains reliable? 



The companion paper is this volume by Kerr addresses some of these questions 12111 . 

3 Weak solutions of Navier-Stokes and a range of scales 

Weak solutions are natural for the global enstrophy 1 1 CO 1 1 ^ because of the properties of projection operators. The 
original argument used by Leray 1 1311 gives us the textbook result from his energy inequality [|4L |5(] . In terms of D\ 
this is 

(D\) T <cRe 3 + 0(T- 1 ) (20) 
where the time average up time T given by (-} r is defined by 



{F(-)) T = -\imsup F(z)dz. (21) 

1 F JO 

To obtain similar results for ||ft)||2„, for m > 1 looks difficult not only because the properties of projection opera- 
tors do not naturally extend to the higher spaces but also because ||ft)||? m does not appear naturally in an energy 
inequality. However, these problems have been circumvented in 11411 . the main result from which will be stated 
below and its very short proof repeated for the benefit of the reader : 

Theorem 1 For 1 < m < °° weak solutions obey 

{D^j^cRei + OiT- 1 ) , (22) 
where c is a uniform constant. 



Proof: The proof is based on a result of Foias, Guillope and Temam [22j] (their Theorem 3.1) for weak solutions. 
Doering and Foias B23I1 used the square of the averaged velocity Uq = L~ 3 (H^HDr t° define the Reynolds number 
Re = UqLv~ 1 which enables us to convert estimates in Gr to estimates in Re. Thus the result of Foias, Guillope 
and Temam [22] in terms of Re becomes 

(h^^ <c N L- l v^Re 3 + 0(T- 1 ) , (23) 
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where 

H N = f r \V N u\ 2 dV = JjV N - l co\ 2 dV , (24) 
and where H\ = fy, \ Vu \ 2 dV = J r \ CO \ 2 d V . Then an interpolation between 1 1 CO \ \ % m and 1 1 CO \ \ 2 is written as 

\\co\\ 2m < c N , m \\V N - l C0\\ a 2 a = , (25) 

for N > 3. 1 1 CO 1 1 2/« is raised to the power A m , which is to be determined. 

■o\\t) T < ^(lIV^fflllf-Hffllli 1 -^ 



, x \aA m {2N-\) [ 
Am I I H TN=T \ 



N,m \ I N / 11 1 



r 

1. 



/ , v ^A m (2W-l) / (l- a) A„, \ 1-I»A„,(2W-1) 

< 4: m \ H N =l ) T (hI^-^ (26) 

An explicit upper bound in terms of /?e is available only if the exponent of H\ within the average is unity ; that is 

{l - a)Am -1 =* A m = ^ = a m { 21) 



2-aAJlN-l) Am -3 

as desired. Using the estimate in (l23l . and d20l i for (Hi), the result follows, c^.m can be minimized by choosing 
N = 3. C3,„, does not blow up even when m = °° ; thus we take the largest value of cf" n and call this c. ■ 

Following the statement of TheoremQ]and motivated by the definition of the Kolmogorov length for m = 1, a 
continuum of length scales X m can be defined thus : 

(D m ) T := (LX m l f am , (28) 

in which case d28l becomes 

LX m l <cRe 3/2a - . (29) 

When m = \, CX\= 2, and thus LAf 1 < cRe 3 / 4 , which is consistent with Kolmogorov's statistical theory [22,23]. 
However, the bounds on X~ 1 become increasingly large with increasing m reflecting how the L 2 '"-norms can detect 
finer scale motions. 

4 A regularity criterion based on the relative sizes of D m and D m+ \ 

Consider two m-dependent constants ci .„, and C2. m and two frequencies GS\ m and GJ2,,,, defined by 

®l, m = (DoOmCi l, ®2,m = ®0<XmC2, m ■ (30) 



In Il2ll , using a standard contradiction strategy on a finite interval of existence and uniquness [0, T*), it was shown 
that for the decaying Navier-Stokes equations the D„, obey the following theorem in which the dot represents 
differentiation with respect to time : 

Theorem 2 For 1 < m < °° on [0, T*) the D m (t) satisfy the set of inequalities 

^Dil-m^^j^y +w 2 , m |, (3i) 

where p m = im(4m +1). In the forced case there is an additive term GJ3 m Re 2 D m . 
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The obvious conclusion is that solutions come under control pointwise in t provided 

D m+l {t)>c Pm D m (t) (32) 

where 

Cp m = [ci, m C2,m] 1/Pm (33) 



This is a numerically testable criterion : if ( l32l holds then the D m must decay in time. It was shown in IU2I1 that 
this can be weakened to a time integral result : 

Theorem 3 For any value ofl< m < °° if the integral condition is satisfied 

'in (1±ZH ) dx > 0, Z m = D m+l /D,„ (34) 

V C 4.m / 

with c 4m = [2 pm ~ l (1 +c\. m c 2 , m )Y'" , then D m (t) < D m (0) on the interval [0, t]. 

Given the nature of c 4m \ 2 it is clear that there must be enough regions of the time axis where D m+ \ > 
(c4, m — \ )D m to make the integral positive. The D„, are easily computable from Navier-Stokes data. Therefore, an 
interesting numerical experiment would be to test : 

1. Whether the D„, are ordered as time evolves such that D„, > D m+ \ or D,„ < D m+ \1 

2. Do the D m cross over from one regime to the other? 

3. How significant are the initial conditions and the Reynolds number in this behaviour? 

5 Body-forced Navier-Stokes equations 

5.1 A critical lower bound on Jq D m dx in terms of Re 



In 1121 12611 the body-forced Navier-Stokes equations were considered in terms of the Grashof number Gr. It is 
more useful to to consider this in terms of the Reynolds number Re. The inclusion of the forcing in <ED modified 
this but requires the introduction of a third frequency GJ3 ,„ = fl7oO£ m C3 .,„ 



t>m < D 3 m I -m,m l^ 1 ) + ®2,m | + ®l,mRe 2 D m , (35) 

where p m — ^m(4m + 1) and y m = \(X m+ \ (m 2 — l) . Let A„, be defined by (2 < A,„ < 6) 

A„, = 3 {8,„(2 + p m y m ) - p m y,„} (36) 

The following result shows that if a critical lower bound is set on / D m d T then D„, will decay exponentially. Note 
that the case m = 1 is excluded : 

Theorem 4 If there exists a value of m lying in the range 1 < m < °°, with initial data [D m (0)] 2 < C„,Re A '", for 
which the integral lies on or above the critical value 

c m (tRe 35 "' + 772) < J D m dx (37) 



pm 



and 8 m and T72 lie in the ranges 
2/3 +p m y m 



<8 m <l, and 772 > T]!^ 5 '"- 1 ' , (38) 

^ ~r Pmym 

then D m {t) decays exponentially on [0, ?]. 

2 For the forced case the definition of £2,„ requires an additive UJq term to act as a lower bound (T2, 26]. 
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Remark: 8„, \ 1/2 for large m so enough slack lies between the upper and lower bounds on / D,„ d x. 
Proof: To proceed, divide by D^ m to write d35l l as 



1 d 
2dt 



D- 2 ) >X m (D-*)-Uh 



-2\ 



X m — GTi n 



D 



m+ 1 



D„ 



(39) 



A lower bound for pQX m dx can be estimated thus : 



/ An+i 
Jo 



dx 



D 



m+1 



d„ 



D 1 



Pm Pm- 2 

D m Pm dx 



< 



Pi" \ Pm 

0\D m J °- d V [J Q D » 



Dm+\ \ 



Pm~ 2 



1/A» 



and so 



./o 



,(['D m+1 dx) pm . 

x m dx>m m t K m+l ' ^ 



(JoO»</t) 



^ -GJ 3 , m ffle z 



(40) 



(41) 



Recall that p m = |m(4m + 1) and 7,„ = jO^+i (m 2 - l) . It is not difficult to prove that < £l™ +l l £li for 
m > 1, from which, after manipulation into the D m becomes 



n <r n^m/^Ymm 2 n a„,/2m 2 
u m ^ u m +l u \ 



and therefore a Holder inequality gives 



f {) D m+ \dx ( j' Q D m dx 



Jo O m rfr " \J* DidT 



(42) 



(43) 



/' 

Jo 



X m dx > QJi m f~ 



(JoD m dx) 



PmYm+2 



(SoDxdx) 



Pm7m 



Inequality d39l integrates to 



[D m {t)f < 



enp {-2 j' Q X m dx} 



[D m (0)} 2 -2(D2. m Jo^p{-2j Q T X m dx'} dx ' 
( |4"T1 > can be re-written as 

(J^t)"'"^ 2 



Jo - ' (^D l dx) p " ,rm 

> c,„t (® Lm Re A "> - m, m Re 2 



(44) 



(45) 



(46) 

having used the assumed lower bound in the theorem and the upper bound of /q-Di dx. Moreover, to have the 
dissipation greater than forcing requires A,„ > 2 so 8,„ must lie in the range as in d38l because 2 < A m < 6. For 
large Re the negative Re 2 -term in d45l > is dropped so the integral in the denominator of d45l > is estimated as 



exp(-2y X m dx' j dx < [2c m aS\. m \ 1 Re A '"(l-exp 



2tBi„c m tRe 



and so the denominator of (|43T > satisfies 

Denominator > [D„,(0)]~ 2 — c?. m ci im (2c m ) _1 Re~ A "' ( 1 — exp 



-2Wi, m c m tRe 



(47) 



(48) 



This can never go negative if [D m (0)] 2 > c\ jn C2, m (2c m ]T l Re A/ ", which means D m (0) < C m Re^ 1 



D m {t) 



\ 




to h t 2 

Figure 1 : A cartoon of D m (t) versus / illustrating the phases of intermittency. The range of 8,„ lies in Theorem|4] The 
vertical arrows depict the region where there is the potential for needle-like singular behaviour & thus a break-down of 

regularity. 

5.2 A relaxation oscillator mechanism for intermittency 

Experimentally, signals go through cycles of growth and collapse Ir27il28il29i[30ll so it is not realistic to expect that 
the critical lower bound imposed in Theorem|4] should hold for all time. Using the average notation (•),, inequality 
d45l l shows that if (D m ), lies above critical then D m (t) collapses exponentially. In Figure 1 the horizontal line at 
Re i5m is drawn as the critical lower bound on (D m ) t . 

Above this critical range, D m (t) will decay exponentially fast. However, because integrals take account of 
history, there will be a delay before (D m ) t decreases below the value above which a zero in the denominator of 
d45l l can be prevented (at fi): at this point all constraints are removed and D m (t) is free to grow rapidly again 
in t\ < / < ?2- If the integral drops below critical then it is in this interval that the occurrence of singular events 
(depicted by vertical arrows) must still formally be considered - if one occurs the solutions fails. Provided a 
solution still exists, growth in D,„ will be such that, after another delay, it will force (D m ) t above critical and 
the system, with a re-set of initial conditions at t%, is free to move through another cycle. Thus it behaves like a 
relaxation oscillator. The vertical arrows in Figure 1 label the region, below critical, where 

J D m dx < c m (f Re 38 '" + 773) . (49) 

It is in this regime where where potentially singular point-wise growth of D m (t) could occur which contributes 
little to the growth of the integral J Q D m dx and which does not drive it past critical. No control mechanism for this 
type of growth is known and so the regularity problem remains open. 

6 Conclusion 

The variables D m , as defined in (0, have proved useful in expressing the Navier-Stokes and Euler regularity 
problems in a natural manner. Their use also poses some interesting questions. For instance, while the i2„, must be 
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ordered because of Holder's inequality, this is not the case with the D m because the a m decrease with m. Theorem 
12 suggests that the regime D m+ \/D,„ > c m , where c m is a constant only just above unity, guarantees the decay of 
D,„ and hence control over Navier-Stokes solutions. In terms of numerical experiments, it would be interesting to 
see, from a variety of initial conditions, which of the two regimes 

D m+ \ > D m or D m+i < D,„ (50) 

are predominant and whether there is a cross-over from one to the other. If so, does this depend heavily on the 
initial conditions, such as the contrasting random or anti-parallel vortex initial conditions? Does it also depend on 
the size of Re? Likewise, do solutions of the Euler equations, when in their intermediate and late growth phases, 
track a singular upper bound as in <TT~8T>? 

Acknowledgement : I would like to thank Darryl Holm of Imperial College London and Bob Kerr of the University 
of Warwick for discussions on this topic. 
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